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(SI . ABSTRACT 

J> ■ A large scale SPH+N-body simulation (GADGET) of the concordance ACDM universe 

£NJ ' is used to investigate orientation and angular momentum of galaxy clusters at z = in 

lO ■ connection with their recent accretion histories. The basic cluster sample comprises the 

3000 most massive friends-of-friends halos found in the 500/i -1 Mpc simulation box. 
Two disjoint sub-samples are constructed, using the mass ratio of the two most massive 
progenitors at z = 0.5 7712/mi (m\ > m^), namely a recent major merger sample and a 
steady accretion mode sample. The mass of clusters in the merger sample is on average 
~ 43% larger than the mass of the two progenitors (mi + 7712), whereas in the steady 
^ ■ accretion mode sample a smaller increase of ~ 25% is found. The separation vector 

Oh| connecting the two most massive progenitor halos at z = 0.5 is strongly correlated with 

the orientation of the cluster at z = 0. The angular momentum of the clusters in the 
\Z{ ' recent major merger sample tends to be parallel to orbital angular momentum of the 

^ \ two progenitors, whereas the angular momentum of the steady accretion mode sample 

is mainly determined by the angular momentum of the most massive progenitor. The 
long range correlations for the major and the minor principal axes of cluster pairs 
extend to distances of ~ 100 ft, _1 Mpc. Weak angular momentum correlations are found 
for distances < 20/i -1 Mpc. Within these ranges the major axes tend to be aligned 
with the connecting line of the cluster pairs whereas minor axes and angular momenta 
tend to be perpendicular to this line. A separate analysis of the two sub-samples 
reveals that the long range correlations are independent of the mass accretion mode. 
Thus orientation and angular momentum of galaxy clusters is mainly determined by 
the accretion along the filaments independently of the particular accretion mode. 

Key words: methods: numerical - methods: statistical - galaxies: clusters: general 
- large-scale structure of the Universe 
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1 INTRODUCTION 



Using a sample of 44 Abell clusters Binggeli (1982) s howed 
that c lusters of galaxies are highly eccentric Trd^Millari et afl 
( 1989) measured the X-ray contours of 49 Abell clus- 
ters and found that most clust ers are clearly flatt ened. 
By means of N-body simulations IWar ren et a ll <ll992l) (see 
also lDubinski fc Carlber Jl99ltlCoTe fc Lacevll99rl demon- 
strated that dark matter halos have triaxial shapes and tend 
to be prolate. Furthermore they showed that the prevailing 
sha pe is supported by the anisotropic velocity dispersion (see 
also lTormenl (^97)). Thus the shape of the halo mirrors its 
dynamical state and is likely to be a long-term feature if no 
major disruptions take place. Subsequently the orientation 



of clusters will be identified with the major axis of its mass 
distribution. 



The angular momentum of galaxy clusters is easily ac- 
cessible in numerical simulations. The orientation of the an- 
gular momentum may record information about the mass ac- 
cretion history and/or tidal interactio ns with the surround - 
ing larg e scale struc ture. Accordin g to lDoroshkevichl 111973th 
iPeebled (Il98dh and lWhitd (I1984T) the primary angular mo- 
mentum of bound objects is due to tidal interaction between 
the elongated proto-structures after decoupling from cosmic 
expansion and before turn-around. More recent studies find 
that the angular momentum of dark matter halos is later 
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modified by the merging hist ory of their building blocks, see 
IVitvitska etafl fcOOfl and lPorciani et alJ i2002ajbh . 

The ali gnment of galax y cluster orientations is contro- 
versial since Binggelj II1982T) published his inaugural study. 
He reported that galaxy clusters are oriented relative to 
neighbours at least for separations smaller than 15 ft~ 1 Mpc. 
Moreover, he found anisotropics in the cluster distribution 
on scales up to 50/i _1 Mpc. The first investigat ions, us- 
ing X -ray contours, led to no llUlmer et alJ ll98Sh or only 
weak jRhee fc Latourll99lll significance of alignment effects. 
HoweveT nChamb^r^et^dTll2002() exploring X-ray position 
angles of 45 clusters found significant alignment of clusters' 
orientations with the connecting line to the next neighbour 
for separations of < 20 ft~ 1 Mpc. Analys i ng a l arge set of op- 
tical data of 637 Abell clusters iPlionij lll994l) found highly 
significant alignment effects on scales below 10 /i _1 Mpc that 
become weaker but extend up to 150 /i -1 Mpc. 

In the literature two slightly different scenarios are pro- 
pose d to explain the alig nment of the orientation of clus- 
ters. |Binne^^.Silk (1979) suggested that tidal interactions 
of evolving proto-cluster systems may lead to the growth 
of anisotropies of clusters and to relative orientation ef- 
fects. However. IWestl il994T) found that clusters grow by ac- 
cretion and merging of surrounding matter that falls into 
the deep cluster potential wells along sheet-like and fila- 
mentary higlT_d«ish^y_j|egtons^_JJsing numerical CDM sim- 
ulations Ivan Haarlem et all (Il997t> demonstrate that clus- 
ters are elongated along the incoming dire c tion of the 
last major merger. Further IWest fc Blakesled i2000T) show 
that Virgo's brightest elliptical galaxies have a remark- 
ably collinear arrangement in three dimensions, which is 
aligned with the filamentary structure connecting Virgo 
and the rich cluster Abell 1 367. In a statistica l investi- 
gation of 303 Abell clusters iPlionis et alJ i2003fi confirm 
the alignments of galaxy members with their parent clus- 
ter orientation as well as with the large-scale environment 
within which the clusters are embedded. Therefore, the clus- 
ter formation is tightly connected with the super-cluster 
network that characterises the large-scale matter distri- 
bution in the universe. High-resolution simulations show- 
ing this effect were described by the V i rgo c ollaboration, 
IColberg et all feOOd) . lOnuora fc Thomasl fcOOCh found a sig- 
nificant align ment signal up to scales o f 30 /i~ 1 Mpc for a 
ACDM model. iFaltenbacher et alJ J20021) and lHopkins et alJ 
( 2005) using large scale dissipationless N-body simulations 
of the ACDM universe report alignment of the orientations 
of galaxy clusters for separations < 100/i _1 Mpc. Improved 
statistics due to a muc h larger simulation volume enable 
iKasun fc Evrard l (|2004l) to find cor relations ranging up to 
200fc~ 1 Mpc. lBasilakos et al.l fcOOrfl using the same simula- 
tion than we do, but an eight times increased mass resolu- 
tion, found alignment between cluster sized halos and their 
hosting super-clusters. 

This work aims to link the vector quantities of galaxy 
clusters (orientation and angular momentum) with their re- 
cent mass accretion history. The basic cluster sample com- 
prises the 3000 most massive CDM halos found by a friends- 
of-friends (FoF) approach. Further the basic sample is di- 
vided into two sub-samples, namely a recent major merger 
sample and a steady accretion mode sample. The subdivi- 
sion criterion is based on the mass fraction of the two most 
massive progenitors at z — 0.5. The masses and positions 



of these progenitors hold information about the pattern of 
infall onto the cluster. We find a strong correlations between 
the infall pattern and the cluster vector quantities. Finally, 
since the infall happens along the large scale filaments we 
analyse long range correlations of these vector quantities be- 
twe en cluster pairs by means of mark correlation fu nctions 
( see IFaltenbacher et alJl2002t IKasun fc Evrardll2004) . 

The paper is organised as follows, ^describes the sim- 
ulation and explains the construction of cluster sample and 
sub-samples. In the impact of different mass accretion 
modes (merger/steady accretion) on shape and angular mo- 
mentum of the clusters is investigated . In S0|the range of 
the alignment of the major, the minor, and the angular mo- 
mentum axes of cluster pairs is measured. [J5j summaries the 
results. 



2 SIMULATION AND CLUSTER SAMPLE 

In subsection H2.ll some information about the simulation 
is provided and H2. 21 describes the construction of the basic 
cluster sample. In the last paragraph ( H2.31 the subdivision 
of the basic cluster sample into a merger and a steady ac- 
cretion mode sample is explained. 



2.1 Simulation 

The simulation is performed using th e public available 
GADGET code (see lSpringel et ahlExill) . According to the 
presently favoured cosmological model, we assume a flat uni- 
verse with the parameters Q m = 0.3, Ha = 0.7 and a Hub- 
ble constant of /iioo = 0.7. The linear power spectrum of 
the density fluctuations has been normalised to erg = 0.9. 
The simulation is realized within a cube of 500/i _1 Mpc 
edge length. The cold dark matter and the gaseous com- 
ponent are represented by 256 3 particles each. The mean 
baryonic mass fraction of the universe is assumed to be 
S7f,/fi m = 0.13, consequently the mass of a single CDM par- 
ticle is rridark = 5.4 x 10 11 /i _1 Mg and the mass of a single 
gas particle is m gaa = 8.1 x 10 10 /i _1 Mq. The gas is sim- 
ulated using smoothed particle hydrodynamics (SPH) with 
no baryonic cooling or heating processes taken into account. 



2.2 Basic cluster sample 

For the identification of cluster halos w e apply a standar d 
friends-of-friends (FoF) algorithm (e.g.. iDavis et af]ll985l) . 
The CDM and the gaseous particle distributions are treated 
separately. The basic cluster is based on the 3000 most mas- 
sive CDM halos found with a linking length o f b = 0.17 x Z, 
wher e I is the mean particle distance (see e.g. Ijenkins et alJ 
l200lf) . The most massive cluster in the basic sample com- 
prises 4201 particles which corresponds to a CDM mass of 
2.2 x 10 15 /i _1 Mq. The smallest object in the sample has 
a CDM mass of 1.1 x 10 14 /i _1 M Q containing 209 particles. 
Subsequently the gaseous halos, identified by an analogous 
FoF approach (b = 0.17 x I) applied to the gas particle dis- 
tribution, are associated with their CDM clusters. We apply 
a linking length of b = 0.17 x I also for redshifts 
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2.3 Mass accretion modes 

ISuwa et alJ J2003f) (see also lRichstone et alll992T) show that 
the formation rate of galaxy clusters in the ACDM universe 
peaks at z ~ 0.4. In the following we want to subdivide the 
basic cluster sample according to the accretion history into 
a recent major merger sample and a steady accretion mode 
sample. Therefore we locate for every cluster the two most 
massive progenitors at z = 0.5, just a short time interval 
before the accretion rate is maximal. These two halos at 
z — 0.5 are identified as being the most massive progenitors, 
which contain the highest numbers of particles also found in 
the current (z — 0) cluster. 

The basic cluster sample is subdivided based on the 
mass fraction of the two most massive CDM progenitors 
mi/m\, where mi and mi denote the most and the second 
most massive progenitor, respectively. The 3000 clusters are 
separated into three disjoint sub-samples, each of them com- 
prising 1000 clusters. The recent major merger sample con- 
tains the clusters with the highest values (0.2 < m^/mi < 1) 
for the progenitor mass fractions, meaning that the two 
progenitors have comparable masses. The steady accretion 
mode sample contains the 1000 clusters with the lowest val- 
ues (< 0.1) for the progenitor mass fractions, meaning that 
the most massive progenitor is by far more massive than 
the second most massive progenitor. Consequently all addi- 
tionally accreted mass must approach within even smaller 
lumps. The 1000 clusters with intermediate progenitor mass 
ratios are omitted from any further investigation, which fo- 
cuses on the impact of the particular accretion mode. 

The upper panel of Fig. Q shows the distribution of 
mass fractions ((mi + mjj/mo) versus mj/mi with the 
steady accretion mode clusters appearing on the left and 
the recent major merger clusters appearing on the right, mo 
denotes the mass of the emerged cluster at z — 0. The inter- 
mediate accretion mode clusters, which are not investigated 
in the following, are located in the middle of this panel (grey 
area). The mean values of the fractions (mi +7712) /mo indi- 
cated by the horizontal lines in the upper panel of Fig.Qare 
~ 24% lower in the merger sample compared to the steady 
accretion sample, whereas mo is ~ 18% higher in the merger 
sample compared to the steady accretion sample. The dif- 
ference in the relative decrease of (mi + mi) /mo and the 
increase of mo indicates a sightly higher additional accretion 
activity in the merger sample. The increase of mo indicates 
that equal mass mergers tend to result in more massive ha- 
los. But they also show slightly enhanced accretion activity 
as indicated by the larger relative decrease of the mean val- 
ues of (mi + mi) /mo compared to the relative increase of 
m . 



3 IMPACT OF ACCRETION MODES ON 
VECTOR QUANTITIES OF CLUSTERS 

In this section we investigate the impact of the pattern of in- 
fall on the principal axes and the angular momentum of cur- 
rent clusters. We show the regularity of the infall of the two 
most massive progenitors in the merger sample and demon- 
strate that an approximate picture of the infall pattern can 
be derived from the information provided by the current 
cluster position and the two most massive progenitor po- 
sitions (at z = 0.5). Subsequently the correlation of these 
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Figure 1. The masses of the current clusters are denoted by mo 
and the masses of the two most massive progenitors at 2 = 0.5 
are denoted with mi and mi, respectively (mi > mi). The up- 
per panel shows the ratio of the combined mass mi + mi divided 
by mo in dependence of mi /mi. In the middle panel the ratios 
of mi /mo and mi /mo (black and gray dots, respectively) are 
shown. The lower panel represents the masses of the current clus- 
ters scaled by the mass of the largest object ~ 2.2 X 10 15 h~ 1 Mq 
at that epoch. The shaded area entitled 'omitted' contains the 
1000 clusters with the intermediate ratios of mi/m\, whereas 
'merger' and 'steady' indicate the sub-samples with the highest 
and lowest mi /mi respectively. The horizontal lines in the up- 
per and lower panels indicate the mean values of the according 
sub-samples, 0.75, 0.64, 0.57 for the upper and 0.09, 0.10, 0.11 
in for lower panel. The dashed line in the lower panel marks the 
boundary of completeness, below this line the second most mas- 
sive progenitors comprise less than 5 particles. 

positions with the shape and the angular momentum of cur- 
rent clusters is examined. 



3.1 Infall pattern 

The following illustration of the infall pattern is based on 
the recent major merger sample (see § I2.3B . The vectors, 
which will be later correlated with the shape and the angular 
momentum of the emerged cluster, are combinations of the 
position vectors of the two most massive progenitors ri and 
ri at 2 = 0.5 and of the cluster's position ro at z = 0. The 
2000 progenitors span a mass range of 1.1 x 10 13 < M < 
4.8 x 10 14 ft _1 MQ which corresponds to a range in particle 




Figure 2. The sketch depicts the vectors used to describe the 
geometry of the accretion process, ri and T2 are the locations 
of the two most massive progenitors at z = 0.5, the vector 
is the displacement vector between theses two and a marks the 
orientation of the cluster at z = (double arrows indicate that 
only the direction not the sign is considered here). The vectors vi 
and V2 are the velocities of the progenitors (compare to Equ. I5Y 

numbers of 24 < N p < 1062. Subsequently we focus on the 
difference vector r<j and the cross product of the velocities 
v x as defined in the following equations. 

r d = ri - r 2 (1) 

V x = Vi X v 2 , (2) 

where 



Vd is the difference vector of the locations of the two most 
massive progenitors. The directions of the mean velocities 
vi and V2 are calculated by subtracting the position of the 
current cluster by the corresponding positions of the pro- 
genitors at z = 0.5. Thus v x is perpendicular to the orbital 
plane of the merging progenitors. In Fig. [21 I'd and v x are 
sketched along with a schematized merging event. 

We only examine the tracks within the recent major 
merger sample to justify the description of the infall pattern 
by ra and v x . Tracking the halos through 33 intermediate 
snapshots in the steady accretion mode is unfeasible due to 
the involvement of low mass (or sparse number of particle) 
halos. Fig. |3] aims to explain the projection and the adjust- 
ment procedures which need to be applied before all the 
tracks can be stacked up. Starting at z = 0.5 we trace the 
2000 most massive progenitors through 33 simulation out- 
puts between z = 0.5 and z — 0. By definition all of these 
tracks have to merge at some point with another track. The 
descendant of an actual halo in the subsequent snapshot is 
in the first instance associated with the halo that contains 
the most overlapping particles. At the same time we check 
whether the masses of the descendant and the actual halo 
are in agreement. For strong deviations (factor of > 10) we 
choose, if present, an alternate descendant showing slightly 
fewer overlapping particles but better mass agreement. 

Fig. [3] shows the projection of two different merging 
events onto the merging plane perpendicular to v x (see 
Eq.|5J. Here two different representations of the same merg- 
ing event are depicted in a row. Each dot represents the 
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Figure 3. Two examples of the adjustment procedure for the 
merging tracks. A row represents the same tracks differently 
stretched. The dots indicate the comoving coordinates of the ha- 
los found in the 33 subsequent snapshots between z = 0.5 and 
z = 0. The tracks are projected onto the plane perpendicular to 
v x (sec Eq. [2J . The tracks on the left panels are shifted to x = 
for z = 0.5 and equal absolute values of the y-coordinates. The 
dashed line connects the origin and the final position of the cur- 
rent cluster. In the panels on the the right a correction for the 
motion in y-direction is applied (i.e. shearing the coordinate sys- 
tem such that the dashed line in the left panels coincides with 
the x-axis) and the merger point of the two progenitor tracks is 
shifted towards the origin. These adjustments do not affect the 
orientation of and v x (see Eqs. ITlandl2t. 

location of the particular halo at a given snapshot. The co- 
ordinate systems in the right panels are rotated in such a 
way that the positions of the two progenitors are located 
on the y-axis (x = 0) with equal separation from the ori- 
gin. The y-coordinate of the most massive progenitor is dis- 
played at positive y-values, consecutively the second mas- 
sive progenitor has negative y-values. In the right panels of 
Fig. [3J the image is stretched in such a manner, that the 
initial coordinates of the progenitors remain unchanged but 
the final cluster position is projected to the x-axis (i.e. the 
dashed line in the left panels coincides with the x-axis in 
the right panels). Finally the merger point is shifted to the 
origin. This adjustments neither change the direction of r<j 
nor the direction of v x . After applying these adjustments to 
all the members of the merging sample the resulting tracks 
are stacked up. Fig. [I] presents face and edge on sight of all 
the merging events. 

As the tracks of particular merging events (Fig.|3J show 
the two progenitors approach each other until they merge in 
a jump-like manner. This behaviour is due to the application 
of the FoF algorithm for the identification of the halos. The 
jumps of ~ 1 h~ Mpc to the merger point reflect roughly 
the radii of the progenitor halos. The separation vector con- 
necting the last distinguishable positions of the progenitors 
is in both cases (upper and lower panels) remarkable aligned 
with rd- This feature can also be seen in the stacked up plot 
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Figure 4. The traces of two most massive progenitors ending in 
the clusters of the merger sample are overlaid after an adjustment 
procedure like sketched in Fig.JJJ has been applied. 



(Fig. The two branches depicting the incoming progeni- 
tor halos show enhanced densities slightly left of x = 0. The 
line joining this two areas is parallel to the y-axis. Thus the 
Vd gives a good approximation of the main infall direction. 
As order of magnitude estimates show the mutual gravi- 
tational attraction of two clusters in the considered mass 
and distance range dominates over the impact of large scale 
gravitational field. Therefore, the two halos are expected to 
approach each other in a fairly radial way, along the line 
joining them. 

The tendency of the encounters to be parabolic can be 
inferred from the average shape of the infall pattern. The 
flatness of the distribution of tracks in the edge on repre- 
sentation (lower panel of Fig.[IJ indicates that v x describes 
the orbital plane in a rather satisfactory manner. These out- 
comes confirm that despite the apparent simplifications the 
vectors and v x are capable of determining the direction 
of infall and the orbital plane, respectively. 
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Table 1. Results of a Kolmogorov-Smirnov test comparing the 
angle distributions indicated in the first column derived from the 
recent major merger sample (second column) and the steady ac- 
cretion mode sample with the random distributions as imaged in 
Fig-OQ 



3.1.1 Impact of infall pattern on clusters' shape 

In this paragraph we investigate the correlations between 
vectors r<j and v x and the shapes of the current clusters. 
The shape of a cluster is determined by the eigenvectors of 
the second moment of mass tensor 9 as given by the following 
equation. 



(3) 



Here mu specifies the particle masses. The position vector 
r of the N CDM particles belonging to the halo is given in 
respect of the centre of mass and the components of the vec- 
tor belonging to the fcth particle are denoted rti for i = 1,3. 
Here we show only the results based on the dark matter 
component, but the behaviour of the gaseous component is 
very similar. Subsequently the three eigenvectors of 8 are 
denoted a, b and c. The eigenvalues o, b and c are ordered 
by magnitude (a > b > c). Hence a indicates the major, 
b the intermediate and c the minor axis. Sin ce cluster ha- 
los show preferentially prolate shape (see e.g. JWarren et all 
Il99l iFaltenbacher et ai]|2002l) a is identified with the ori- 
entation of the cluster. 

The distributions of the angles between the principal 
axes and iv are shown in Fig. |S] The diagonal represents 
the cumulative distribution of the cosine of the angles be- 
tween two randomly orientated vectors. A curve lying below 
the diagonal indicates an correlation between vectors, i.e. 
the vectors have a tendency to be parallel more often than a 
random sample. A curve lying above the diagonal indicates 
a anti-correlation between the vectors, i.e. the vectors have 
a tendency to be more often perpendicular than in a ran- 
dom sample. The thick solid, dashed and dotted lines depict 
the cumulative distributions of the angles between and 
the principal axes a, b and c as indicated. The left panel 
shows the results for the recent major merger sample and 
the right panel those for the steady accretion mode sample. 
The dashed lines that represent the angles between the ma- 
jor principal axis a and rd deviate most significantly from a 
randomly generated distribution (diagonal) . The major axis 
(a) tends to be parallel to r<j. The result of a Kolmogorov- 
Smirnov (KS) test, comparing these distributions with the 
random distribution, is given in Tab. Q The opposite is 
found for the distributions associated with the two smaller 
principal axes. They clearly show a trend to be perpendic- 
ular to Td- This behaviour is expected since the principal 
a, b and c are perpendicular to each other. Therefore, if a 
vector distribution is strongly aligned with one of these axes 
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Figure 5. The cumulative distribution of the absolute values 
of the cosine between the difference vector at z = 0.5 and 
the the principal axes a (dashed), b (dotted) and c (solid) of the 
cluster at z = is shown. The diagonal indicates the shape of 
a random distribution. The left panel gives the result from the 
merger sample and the right panel shows the outcome for the 
steady accretion mode sample. 



Figure 6. The cumulative distribution of the absolute values of 
the cosine between the cross product v x of the two infall veloci- 
ties vi and V2 and the the principal axes a (dashed), b (dotted) 
and c (solid) of the cluster at z = is imaged. The diagonal indi- 
cates the shape of a random distribution. The left panel gives the 
result from the recent major merger sample and the right panel 
shows the outcome for the steady accretion mode sample. 



it has to be perpendicular to the other two principal axes. 
The different deviations of ip(h, r<j) and if>(c,rd) from ran- 
dom expresses a strong confinement of to the a-b plane 
(still the alignment with a is by far the strongest). The KS 
probabilities that the actual distributions are accidental re- 
alizations of a random process fall below 10~ 45 . The gen- 
eral trends match for the distributions related to the two 
different samples (merger/steady). The distributions of the 
merger sample show somewhat enhanced deviations from 
random. The shape of clusters which underwent a recent 
major merging event is obviously dominated by the orbital 
parameters of the two progenitors, expressed by the vector 
r<j. However, since the results of the steady accretion mode 
sample show qualitatively equal behaviour, the orientation 
of clusters must be determined in a more general sense by 
the coherent mass accretion along the filaments of the cos- 
mic network. 

As discussed above, § 13.11 v x is perpendicular to the 
plane of infall and parallel to the orbital angular momen- 
tum of the two progenitors. In Fig. [5] the solid, dotted and 
dashed lines correspond to the distributions of the angles 
between v x and the minor, the intermediate and the major 
axis. Again the deviation from a randomly generated angle 
distribution (diagonal) is obvious. The most significant de- 
viation from random is seen for ip(a, v x ), whereas ip(h, v x ) 
are almost compatible to random orientations. The KS test 
results are displayed in Tab. |5] The shape of the distribu- 
tions indicate, that the major axis is perpendicular to v x , 
whereas the minor axes tend to be in parallel with v x . As 
mentioned above the distributions related to the different 
principal axes are not independent. A vector distribution 
which is strongly aligned with one of these axes is expected 
to be perpendicular to the other two axes. This is apparently 
not the case for i/)(b,v x ). This behaviour can be explained 
if one takes the uncertainty in the determination of the clus- 
ter axes into account, which accumulates in the distribution 
of b. On the one hand, if the cluster is prolate (with eigen- 
values a > b ~ c) then b and c are degenerate, on the other 
hand, if the cluster is oblate (a ~ b > c) then b and a are 
degenerate. In total both effects result in an almost ran- 
dom distribution for ^(b, v x ). The distributions associated 
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Table 2. Results of a Kolmogorov-Smirnov test comparing the 
angle distributions indicated in the first column derived from the 
recent major merger sample (second column) and the steady ac- 
cretion mode sample with the random distributions as plotted in 
FigEl 

with the merger events show a more prominent deviation 
from the random for tp(a., v x ) and ip(c,\ x ) and a slightly 
more random like distribution for ip(h, v x ) compared to the 
steady accretion mode sample. However (as seen before in 
the Td distributions) the similarity of the v x distributions 
indicates that the mass accretion along the filaments is the 
predominant factor in the determination of the orientation, 
no matter which accretion mode (merger/steady) is consid- 
ered. 

3.1.2 Impact of infall pattern on clusters' angular 
momentum 

The angular momentum is closely related to the dynam- 
ics of the halo. However, since the spin parameter A of ha- 
los is roughly 0.05 the ord ered rotational support against 
gravity is only abo ut 5% dSteinmetz fc Bartelmannl Il995l : 
ICole fc Lacevlll996l) . Thus the correlation analysis using the 
direction of the angular momentum is expected to be blurred 
by this effect. The angular momentum for a FoF cluster is 
computed in the following way. 

1 = 2J ™i= r i x v ' (4) 

The r denotes the particle positions relative to centre of 
mass, rrifc is the particle mass and N is the number of par- 
ticles belonging to the halo. As discussed above the vector 
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Figure 7. The cumulative distribution of the absolute values 
of the cosine between; 1. the angular momentum vector lo at 
2 = and the cross product v x of the two infall velocities vi 
and V2 at z = 0.5 (solid line), 2. lo and the minor principal axis 
c (dashed-dotted), 3. lo and the angular momentum vector li of 
the most massive progenitor at z = 0.5 (dotted line), 4. li and 
the the angular momentum vector I2 of the second most massive 
progenitor at z = 0.5 (dashed line). The diagonal indicates the 
shape of a random distribution. The left panel gives the result 
from the recent major merger sample and the right panel shows 
the outcome for the steady accretion mode sample. 
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< 10" 45 
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-1 


7.5 x 10" 1 
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1.0 x 10" 


35 


2.3 x 10~ 17 



Table 3. Results of a Kolmogorov-Smirnov test comparing the 
angle distributions indicated in the first column derived from the 
recent major merger sample (second column) and the steady ac- 
cretion mode sample with the random distributions as presented 
in Fig|7| 

v x is perpendicularly orientated to the infall plane. If the 
angular momentum lo of the emerged cluster at z — is 
dominated by the orbital angular momentum of two most 
massive progenitors, then a correlation between v x and lo 
is expected. Therefore we analyse subsequently the distri- 
bution of the angles between lo and v x . Additionally we 
consider the angular momenta of the two progenitor halos 
li and I2, where li corresponds to the angular momentum 
of the most massive and I2 to the second most massive pro- 
genitor. In Fig. [T] the various angle distributions connected 
to the different angular momenta are displayed. Since we 
have found a strong correlation between v x and c (minor 
axis of the cluster) we also show ^(lo, c) using dashed-dotted 
lines. On the left hand side the results for the merger sam- 
ple are depicted and the outcome for the steady accretion 
mode sample is displayed on the right panel. The accord- 
ing KS test results comparing these distributions with the 
random distributions (diagonals in Fig. [7J are presented in 
Tab. |3 The distributions between the angular momenta of 
the two progenitors ^(li, I2) agree with random orientations 
for both samples. V(lo,c) also shows similar behaviour, in- 
dicating that lo tends to be parallel to c. Here the KS prob- 



abilities are enhanced for the merger sample, but the quali- 
tative trend is similar in the steady accretion mode sample. 
The most obvious differences between the two samples are 
seen in ?/>(lo,v x ) and i/j(1o,1i). These distributions show op- 
posite behaviour for the merger and the steady accretion 
mode sample. In the merger sample the KS probabilities for 
i/i(lo,v x ) are very low, indicating that this distribution is 
obviously incompatible with a random distribution. Instead 
tp(lo, li) is almost in agreement with a random distribution. 
The opposite picture arises in the steady accretion mode 
sample. Here ip(lo, li) is clearly not in agreement with a ran- 
dom distribution, whereas ^>(lo,v x ) can't be differentiated 
from random. 

We find the effect of steady accretion on the angular 
mom entum is in line with previous work (i.e. lVitvitska et alJ 
2002). If there is no major merger, the angular momentum of 
the most massive progenitor will be aligned with the angu- 
lar momentum of the descendant cluster. As seen in Fig. Q 
the most massive progenitor in the steady accretion sample 
is orders of magnitude larger than all other infalling halos. 
This explains the result for t/j (lo , li). The small tendency of 
an alignment bet ween lp and v x may have the same origin 
as the finding bv iBensonl j200Fj) . that merging dark matter 
satellites show a weak tendency to have velocities perpen- 
dicularly orientated to the angular momentum of the main 
halo. 

In the merger sample the orbital angular momentum of 
the two most massive progenitors shows a strong impact on 
the angular momentum of the descendant. A decided differ- 
ence between i/>(lo,v x ) and a random distribution can be 
seen. The KS probability of 2.1 x 10~ 13 strongly discounts 
any random process as being the origin of this distribution. 
The tendency for lo and v x to be parallel confirms the pic- 
ture that a large fraction of the descendant's angular mo- 
mentum was gained from the orbital angular momentum of 
the two most massive progenitors. This result is in agree- 
ment with the detail ed study of the angula r momenta of 
dark matter halos bv IVitvitska et all i2002l) . Whereas the 
correlation of li and lo is in agreement with random orien- 
tations as can be seen by the rather high value of the KS 
test. 



4 CORRELATIONS IN ORIENTATION AND 
ANGULAR MOMENTUM 

As shown above orientation and angular momentum of cur- 
rent clusters are tightly connected to the infall pattern. 
The dominant amount of mass is accreted along filaments. 
Therefore we employ the technique o f mark correlation 
functi o ns (see [Stoy an fe Stoyanl il994) : iFaltenbacher et alJ 
J20021) : iKasun fc Evrardl fcooli ) to investigate the differ- 
ence in the large scale correlations between the basic sample 
and the two accret i on mo de sub-samples (merger/steady). 
IFaltenbacher et alJ <|2002T) used the following equation to 
study filamentary alignment. 

^v(r) = i(|v 1 -f| + |v 2 -f|} p (r), (5) 

Here v is a normalised vector quantity (|v| = 1), r is the dis- 
tance vector between a pair of clusters with absolute value 
r. The normalised distance vector is denoted by r = r/r 
and ( ) p denotes the average over all pairs with separation 
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r. Due to the spatial symmetry the absolute values of these 
scalar products are used. J-(r) quantifies the filamentary 
alignment of the vectors vi and V2 with the line connecting 
both clusters. T(r) is proportional to the mean value of the 
cosine of the angles between vi and the vector f connect- 
ing the pairs. Random orientations imply J-(r) — 0.5. Since 
the distribution of angles is not of Gaussian shape around 
a given mean value, it would be wrong to apply Gaussian 
statistics to estimate the the errors. Thus, to investigate 
the la confidence intervals we per f orm a non-parametric 
Monte Carlo test <Besag fc Digglej Il977l : iGottlober et alJ 
l2002l:lKerscher et all200lHKerscherll998ft . Here the null hy- 
pothesis is "no deviation from random angles" . We simulate 
this null hypotheses by keeping the positions of the halos 
fixed and randomly re-assigning the orientations to the ha- 
los. This procedure is repeated 99 times. The rms-deviation 
for these realisations gives a measure for the la confidence 
interval for the given data set. These intervals are indicated 
by the shaded areas around the expected values for random 
distributions in Fig. |S| Thus, if J-(r) lies outside of this in- 
trinsic scatter region around 0.5 the deviation from random 
orientations is confirmed at a la confidence level. In addi- 
tion to the alignment of the CDM compon ent we present the 
alignm ent signa ls for the gaseo us halos. iKazantzidis et al] 
J2004I) (see also lDubinsk1ll99^ pointed out that the elon- 
gation of halos is reduced due to radiative cooling and subse- 
quent condensation of baryonic matter at the cluster centre. 
Since we use non radiative SPH simulations this effect is not 
represented. However according to Eq.[3]tlie second moment 
of mass tensor 6 and thus the orientations of the principal 
axes are dominated by the mass distribution on the outskirts 
of the clusters. Moreover, we are mainly concerned with the 
orientation of the halo, not the actual axial ratios, and the 
orientation is not affected by dissipative effects. 

4.1 Correlation of clusters' orientations 

We have shown in that the different accretion modes do 
not result in a qualitatively different impact on the shape 
of current cluster. However the correlation between shapes 
and infall pattern are reduced in the steady accretion mode 
sample. Here we separately investigate the filamentary cor- 
relation of the orientation of cluster for the different sam- 
ples (basic, merger, steady) to see whether this quantitative 
differences influence the signal for large scale correlations. 
Replacing v in Eq.|S]with the major (and subsequently the 
minor principal) axes results in 

^p(r) = i<|pi-f| + |p 2 -f|) p (r) • (6) 

Here pi and p2 denote the principal axes (major or minor) of 
a cluster pair at z — and the associated separation vector 
r. The analysis is done for the whole cluster sample and for 
the two sub-samples (merger/steady) separately. The results 
for the alignment of the major and the minor axes are shown 
in upper and middle panel of Fig. [H] 

F{r) > 0.5 indicates that the axes tend to be aligned 
with the connecting line of the clusters. For the basic clus- 
ter sample, comprising the 3000 most massive clusters in 
the simulation volume, this signal for the major axes ex- 
tends up to 1 00/i~ 1 Mpc (left panel of Fig. |S| compare also 
to Fig. 3b in iFaltenbacher et alJ 12002). The signal for the 



merger and the steady accretion mode samples extends up 
to 80/i _1 Mpc. There is no distinctive difference between the 
results for the two sub-samples. The somewhat shorter range 
seen in both sub-samples appears to be a result of the lower 
statistics (only 1000 members per sub-sample) rather than 
caused by the different accretion modes. 

3~(r) < 0.5 indicates that the axes tend to be perpen- 
dicular to the connecting line of the cluster pairs. In the 
basic cluster sample the signal for the minor axis extends 
up to 100ft~ 1 Mpc (middl e left panel of Fig. Et. Thi s result 
is in good agreement with lBailin fc Steinmetz| i2005l) . where 
a strong tendency for the minor axes of halos to lie perpen- 
dicular to the large scale filaments is found. The CDM com- 
ponents of the recent major merger sample (middle panel 
of Fig. |HJ and the steady accretion mode sample show a 
tendency to be perpendicular to the connecting line for sep- 
arations up to 70/i _1 Mpc. Here again the shorter correla- 
tion ranges for the sub-samples are likely be caused by the 
smaller sample size. 

The comparison between the signals for the CDM com- 
ponent (thick solid line) to the gaseous component (thin 
dotted line) generally show good agreement. However, the 
small differences in the correlation range are likely to mirror 
the different physical pro perties of the two comp onents. Us- 
ing X-ray position angles Ichambers et aljl2002l found that 
clusters with separations d n < 20/i _1 Mpc tend to point to- 
wards their neighbours. Their sample comprised 45 clusters 
thus the statistics are fairly low compared to our analysis. 

Previously we have shown that the infall of halos is 
correlated with the orientation of the cluster and is inde- 
pendent of accretion mode. In this section we demonstrate 
that the orientations of clusters are correlated and indepen- 
dent of accretion mode. Since the steady accretion mode 
sample shows just as much correlation as the merger mode 
sample all progenitors must be infalling along the same rel- 
ative direction. From this, we conclude that the alignment 
is strongly supported by infall along filamentary structure. 
The correlation length (< 100/i _1 Mpc) for cluster align- 
ment that we find in our simulatio n is in agreement with 
the findings of [Hopkins et al] (2005). 

4.2 Alignment of clusters' angular momenta 

In § 13.1.21 we have analysed the effects of different accre- 
tion modes on the angular momentum of current clusters. 
We find a different behaviour in the merger and the steady 
accretion mode samples. To see, whether these differences 
have impact on the large scale correlations, we investigate 
the correlation of the intrinsic angular momenta of cluster 
pairs in the whole sample, as well as in the two sub-samples. 
Therefore the angular momenta is used as vector mark in 
equation J5J. 

*k(r) = i(|Ia-?| + [l*-?|)p(r) (7) 

Here l a and If, denote the normalised angular momenta of a 
cluster pair separated by the spatial vector r. The result of 
this study is displayed in the lower panels of Fig. |H| Similar 
to the result for the correlation of the minor axes (middle 
panels Fig. the angular momenta of the CDM component 
show the trend to be perpendicular to the connecting lines. 
Thus they tend to be perpendicular to the filaments as well. 
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Figure 8. The upper panel shows the mark correlation function, using the clusters' major axes as vector marks. The solid lines display 
the signal for the CDM component. The dotted lines show the behaviour of the gas. The horizontal solid line gives the expected values for 
a random distribution of the orientations. The shaded area is a measure of the intrinsic scatter of the signal. The three panels from right 
to left image the results for the whole, the merger and the steady accretion sample, respectively. In the middle panels the correlation for 
the minor axes of the clusters are shown and the lower panels represent the signal obtained for the angular momentum. 



The signal reaches only out to separations of ~ 20/i~ 1 Mpc 
for the complete sam ple. This result is in quali tative agree- 
ment with Fig. 7b in iFaltenbacher et alJ i2002t) . 

The CDM signal in the recent major merger sample is 
hardly distinguishable from random fluctuations. Further- 
more the deviation from random is far smaller than seen 
in the correlations of the principal axes. As mentioned in 
the beginning of H3.1.2I coherent rotation supports the halo 
against gravitational collapse only at a 5% level, thus a low 
deviation from random may be expected. However, notable 
is the long correlation range of the steady accretion mode 
sample (~ 60 ft _1 Mpc) compared to the almost disappearing 
signal for the merger sample. According to the established 
picture of halos' mass growth they go through a merger 
phase (ra pid mass growth) prio r to a steady accretion phase 
(see e.g. IWechsler et"afl 12002). Therefore, on average the 
merger sample is younger than the steady accretion sample. 
The results from the above combined with those from the 
previous section lead to the conclusion that infall along the 
filaments tends to torque the halo's angular momentum into 
a perpendicular orie ntation with respect to those filaments. 
IColberg et al.l (2005) finds a trend that more massive clus- 
ters have more associated filaments. Fig. indicates a rela- 
tive increase of clusters' masses in the major merger sample. 
Moreover, the mass accretion in addition the two most mas- 
sive progenitors is maximal in the major merger sample. The 
rapid accretion in the merger sample takes place along vari- 



ous filaments which causes violent changes in the the angular 
momentum. In a later state, when the smaller filaments are 
exhausted, the infall direction is more uniform. According to 
this scenario it is expected that the merger sample will set- 
tle down in time and reflect the same sort of correlations as 
the current steady accretion sample. We plan on addressing 
this is sue in a soon to be forthcoming paper lAllgood et alJ 
J2005tl . 

The slightly different behaviour of the gaseous compo- 
nent indicated by the thin dotted lines in the lower panel of 
Fig. |H| is likely to be caused by differences in the underly- 
ing dynamics. The gaseous component mainly deviates for 
cluster pair separations < 20/i~ 1 Mpc from the dark matter 
distribution. 



5 SUMMARY 

We use large scale SPH+N-body simulations within the con- 
cordance ACDM model to measure the mass accretion histo- 
ries and the correlation of orientation and angular momen- 
tum of galaxy clusters. The two most massive progenitors at 
z = 0.5 are located for the 3000 most massive current clus- 
ters. Based on the mass fractions of these two progenitors 
milm\ (mi > rwi) two sub-samples were constructed, a re- 
cent major merger sample containing the 1000 clusters with 
highest m-zlrn\ values and a steady accretion mode sample 
with the 1000 lowest 7712/mi values. 
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The mean values of (mi + 7712) /mo, where mo denotes 
the current mass of the particular cluster, for the recent 
major merger sample and the steady accretion mode sample 
are 0.57 and 0.75, respectively. The clusters which experi- 
enced recent major mergers tend to accrete even more mat- 
ter in addition to mi + ni2 than silently evolving clusters. 
Major mergers predominantly happen to clusters showing 
averagely enhanced mass accretion. Moreover these clus- 
ters tend to have higher masses than their non merging 
counterparts. This o utcome can be related to the work by 
IColberg et al.l (l2005h . They demonstrate that more massive 
clusters tend to have more filaments approaching them. Con- 
sequ ently the i r over -all accretion activity may be enhanced. 
Also iBensonl 120051) reports strongly enhanced merging ac- 
tivities for the largest members in his sample of N-body 
clusters. 

The reconstruction of the halo tracks in the merger 
sample reveals that the infall pattern are not complex. The 
merging takes place along the vector iv joining the two most 
massive progenitor halos at z — 0.5. Despite the relatively 
high additional mass accretion the tracks of the two most 
massive progenitors are reasonably well confined within the 
plane defined by v x (Eq. |5J. The tracks do not twine each 
other, instead the two most massive progenitors approach 
each other parallel to the line joining them z = 0.5. This 
behaviour is expected since the mutual gravitational attrac- 
tion of two clusters in the considered mass and distance 
range dominates over the impact of large scale gravitational 
field. 

We find that the orientation and the angular momen- 
tum of current clusters are tightly correlated with the vec- 
tors Yd and v x • Td tends to be parallel to the major axis and 
perpendicular to the minor axis of the cluster. The opposite 
behaviour is found for v x . Both accretion mode sub-samples 
show similar correlations. This indicates that the current ori- 
entation is a result of mass flows along the filaments onto 
the cluster, independent of accretion mode. It is likely that 
the accretion flow determines the velocity structure of the 
cluster. If this is the case our results agree wi th alignment o f 
posi tion and velocity princip al axes found by Tormen ( 1997) 
and lKasun fc Evrardl 120041) . However, the angular momen- 
tum behaves differently in the two sub-samples. The angular 
momenta lo of the clusters in the merger sample show the 
tendency to be parallel to v x , whereas the orientation of the 
intrinsic angular momenta of the most massive progenitors 
li at z = 0.5 has a minor influence on the angular momen- 
tum at z — 0. The opposite picture arises in the steady ac- 
cretion mode sample. Here lo and li are correlated, whereas 
the angle distribution between lo and v x is in agreement 
with random orientations. 

The major axes of cluster pairs in the basic sample tend 
to be parallel to their connecting line for distances up to 
~ 100/i _1 Mpc. Correlations up to ~ 80/i -1 Mpc are found 
in the steady accretion and the merger sample. The reduced 
extension of the correlation is likely caused by the reduced 
statistic in the sub-samples. The minor axes of cluster pairs 
tend to be perpendicular to their connecting line. The exten- 
sion of these correlations is similar to the findings for the ma- 
jor axes. The appearance of large scale correlations of clus- 
ter or ie ntations has al s o bee n showjiJjyjRd^nbacher^^dJ 
<2002h : iHopkins et all (l2005h and iKasun fc Evrardl 12004) 
The mass accretion along the filaments appears to be re- 



sponsible for the signal, since the accretion mode shows no 
impact on the correlation signal. 

Weak correlations are found for the angular momenta 
in the complete and the steady accretion mode sample. The 
angular momenta tend to be orthogonal to the connecting 
lines of cluster pairs up to ~ 20ft _1 Mpc for the complete 
and up to ~ 60 fo _1 Mpc for the steady accretion mode sam- 
ple. The orthogonality between angular momentum and fil- 
aments is a result of mass accret ion along these filaments 
(see also lBailin fc Steinmet J2005T) . The signal is stronger for 
the steady accretion sample because the halos in this sample 
have already gained major parts of their angular momentum 
from infall along the filaments. The merger sample is still in 
a phase of accretion and it is expected that additional infall 
along the filaments will bring the merger sample into better 
agreement with the steady accretion sample. The correla- 
tions in the orientations and the angular momenta for the 
gas and the dark matter components show very similar be- 
haviour, supporting the search for alignment effects in large 
X-ray samples of clusters of galaxies. 
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